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Scientific  Accomplishments 

This  final  report  summarizes  research  accomplished  during  July  19,  2006  to  July  31, 
2010  on  developing  and  applying  first  principles  quantum  mechanics  methods  to  describe 
mechanical  wear  of  chrome -plated  steel  gun  tubes.  During  this  grant,  we  published  seven 
articles  and  gave  15  invited  talks  and  6  contributed  presentations  (talks  and  posters)  at 
conferences  under  ARO-DRI  sponsorship.  To  accurately  calculate  from  first  principles  the 
mechanical  properties  of  oxides  that  contribute  to  gun  tube  wear  and  erosion,  the  electronic 
structure  of  the  material  must  be  described  properly.  The  first  half  of  this  grant  was  devoted  to 
developing  a  quantum  mechanics  method  that  was  both  accurate  for  such  materials  and  non- 
empirical  so  that  it  could  be  predictive.  This  entailed  development  of  our  ab  initio  DFT+U 
theory,  published  in  Physical  Review  B  and  the  Journal  of  Chemical  Physics.  We  next  published 
an  Invited  Perspective  article  in  Science  on  the  status  and  challenges  in  quantum  mechanics 
based  materials  modeling,  an  article  in  the  Journal  of  Mechanics  and  Physics  of  Solids  that 
discusses  the  response  of  chromia  to  tensile  loads  across  varying  length  scales,  and  an  article  in 
Acta  Materialia,  which  derives  a  method  for  predicting  shear  strengths  across  multiple  length 
scales  and  predicts  the  response  of  chromia  to  shear  forces.  In  2010,  we  published  two  papers 
that  use  the  universal  scaling  laws  we  derived  for  tensile  and  shear  deformations  to  compare  the 
mechanical  response  of  three  iron  oxides,  Fe2C>3,  FeO,  and  Fe3C>4,  across  a  variety  of  length 
scales.  All  of  these  papers  employed  our  ab  initio  DFT+U  theory  that  provides  the  correct 
electronic  structure  for  chromia  and  the  iron  oxides,  producing  accurate  mechanical  properties. 

A.  Ab  Initio  DFT  +U  Method 

Transition  metal  oxides  such  as  Cr9C>3,  FeO,  and  Fe?C>3  are  examples  of  a  larger  class  of 
materials  known  as  strongly  correlated  electron  materials  (SCEMs),  namely  those  for  which  a 
mean-field  description  of  electron  correlation  -  such  as  that  provided  in  density  functional  theory 
(DFT)  -  completely  fails.  In  conventional  DFT  calculations,  a  model  is  chosen  for  the  electron 
exchange-correlation  (XC)  density  functional,  since  the  exact  XC  functional  is  unknown.  By 
construction,  the  model  necessarily  approximates  exchange  interaction  energies.  In  many  cases, 
the  typical  approximations  [the  local  density  approximation  (LDA)  or  the  generalized  gradient 
approximation  (GGA)]  work  extremely  well  for  a  variety  of  ground  state  properties.  However, 
there  exist  classes  of  materials  and  phenomena  for  which  the  LDA  and  GGA  simply  do  not 
contain  the  correct  physics.  Here  we  are  interested  in  one  of  these  cases,  namely  the  case  in 
which  the  material  has  multiple  open-shell  valence  electrons  occupying  tightly  contracted 
orbitals.  The  open  shell  d-electrons  of  the  transition  metal  ions  are  treated  improperly  by  the 
LDA  and  the  GGA,  because  their  approximate  intra-atomic  exchange  energies  do  not  exactly 
cancel  the  spurious  Coulomb  repulsion  of  each  open-shell  d-electron  with  itself,  leading  to  the 
so-called  “self-interaction  error  (SIE)”.  The  consequences  of  the  inexact  cancellation  by 
approximate  DFT  functionals  are  quite  dire:  not  just  quantitative  but  qualitative  failure,  with 
semiconductors  sometimes  predicted  to  be  metals  (e.g.,  NiO  and  FeO).  Why  does  the  electronic 
structure  go  so  horribly  wrong?  To  minimize  this  spurious  net  Coulomb  repulsion,  the  open-shell 
electrons  inappropriately  delocalize,  and  delocalized  electron  distributions  are  metallic.  So  the 
physics  is  ruined  by  the  SIE.  Why  is  the  failure  so  dramatic  just  in  the  cases  enumerated? 
Because  they  each  have  large  numbers  of  large  intra-atomic  exchange  interactions,  each  of  which 
is  in  error  and  the  error  is  cumulative. 
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By  contrast,  this  spurious  self-Coulomb  repulsion  is  exactly  cancelled  by  the  exact 
exchange  energy  in  Hartree-Fock  (HF)  theory.  Thus,  a  solution,  as  first  articulated  by 
Lichtenstein  and  coworkers,1  comes  in  the  form  of  treating  the  intra-atomic  Coulomb  and 
exchange  interactions  properly  with  HF  theory,  while  treating  the  rest  of  the  system  with  DFT. 
This  is  the  essence  of  the  DFT+U  theory,  which  has  been  shown  to  provide  accurate  descriptions 
of  the  ground  state  electronic  structure  of  SCEMs.  Are  there  drawbacks  of  DFT+U  theory? 
There  are  two:  (i)  a  choice  must  be  made  for  the  choice  of  “U”,  actually  U-J,  the  difference 
between  averaged  intra-atomic  Coulomb  (U)  and  exchange  (J)  energies.  The  selection  is  usually 
empirical  (based  on  comparison  with  experimental  observables  on  known  materials),  which  is 
problematic  both  formally  and  practically  -  what  if  one  wants  to  study  a  brand  new  material  for 
which  no  measurements  have  been  made?  Another  approach  involves  constrained  DFT 
calculations  to  determine  U-J,2  but  this  has  the  obvious  drawback  of  establishing  U-J  from  the 
very  theory  one  is  trying  to  correct,  (ii)  Once  chosen,  U-J  is  fixed  for  the  material  at  hand.  While 
this  is  sufficient  in  many  instances,  circumstances  exist  where  a  dynamical  treatment  of  U-J  is 
really  needed.  For  example,  metal-insulator  transitions  are  not  treated  properly  with  DFT+U, 
since  HF  does  not  treat  metals  properly  (generating  spurious  charge  density  waves,  etc.).3  For 
such  phenomena,  use  of  the  dynamical  mean  field  theory  (DMFT)  is  appropriate,  but  such 
phenomena  are  not  of  concern  in  the  present  work. 

Our  contribution  to  DFT+U  theory  addresses  the  first  drawback.  We  developed  a  scheme 
to  determine  U-J  fully  non-empirically  and  without  SIE  from  ab  initio  unrestricted  Hartree-Fock 
(UHF)  calculations  on  electrostatically  embedded  clusters.  We  derived  a  relationship  between  U 
and  J  and  the  Coulomb  and  exchange  integrals  evaluated  using  the  full  set  of  occupied  UHF 
molecular  orbitals  (MOs)  for  the  electrostatically  embedded  clusters.  We  evaluate  U-J  using 
increasingly  larger  clusters  representing  the  crystal  of  interest,  where  the  clusters  are  embedded 
in  very  large  point  charge  arrays  to  represent  the  Madelung  potential  of  the  extended  crystal. 
Cationic  nearest  neighbors  to  the  anionic  atoms  on  the  surface  of  each  cluster  are  not  described 
as  point  charges,  but  instead  are  represented  by  effective  core  potentials,  in  order  to  account  for 
core-valence  electron  Pauli  repulsion  between  the  valence  electrons  of  the  surface  anions  and  the 
neighboring  cation  cores.  The  cluster  size  is  increased  until  we  reach  converged  values  of  U-J. 
The  technique  has  been  used  to  evaluate  ab  initio  values  of  U  and  J  for  Cr203,  FeO,  and  FC2O3 
thus  far.  Unlike  the  empirically-derived  values  for  these  materials,  the  parameters  we  derive 
exhibit  physically  reasonable  periodic  trends.  For  example,  the  value  of  U-J  for  Cr3+  is  smaller 
than  the  corresponding  value  for  Fe3+,  which  makes  sense  since  Fe’s  orbitals  are  more  contracted 
than  Cr’s  due  to  Fe’s  larger  effective  nuclear  charge  (Fe  is  less  well  screened  than  Cr). 
Secondly,  the  value  of  U-J  for  Fe2  is  smaller  than  for  Fe3+,  which  again  makes  sense  since  Fe3+ 
has  more  contracted  orbitals  than  Fe2  .  By  contrast,  empirically  derived  values  for  U-J  give 
exactly  the  opposite  (unphysical)  trends.  DFT+U  calculations  using  the  ab  initio  parameters 
yield  bulk  electronic  and  structural  ground  state  properties  in  good  agreement  with  experiment. 
As  such,  this  approach  provides  the  means  to  perform  accurate  and  fully-predictive  DFT+U 
calculations  of  SCEMs.  Two  papers  have  appeared  on  this  work:  N.  J.  Mosey  and  E.  A.  Carter, 
“Ab  initio  Evaluation  of  Coulomb  and  Exchange  Parameters  for  DFT+U  Calculations,  ”  Phys. 
Rev.  B,  76,  155123  (2007)  and  N.  J.  Mosey,  P.  Liao,  and  E.  A.  Carter,  “Rotationally -Invariant 
ab  initio  Evaluation  of  Coulomb  and  Exchange  Parameters  for  DFT  +  U  Calculations,  ”  J. 
Chem.  Phys.,  129,  014103  (2008). 


2 


B.  Response  of  O2O3  to  Tensile  Loads 

Chrome-plated  gun  tubes  will  have  a  layer  of  chromia  (C^Cb)  on  them  upon  exposure  to  moist 
air.  Understanding  the  mechanical  response  of  chromia  will  give  some  insight  into  wear 
properties  of  the  gun  tube.  Hence,  we  first  evaluated  with  our  ab  initio  DFT+U  theory  various 
mechanical  properties  associated  with  the  fracture  of  chromia  by  applying  tensile  loads  along  the 

[0001]  and  [Oil  (3/2  )(a/c)“  2]  directions.  Loading  along  these  directions  led  to  formation  of 

the  most  stable  surfaces  and  therefore  we  presume  these  directions  to  be  the  most  likely  to  lead 
to  fracture.  Tests  of  this  assumption  provided  justification  for  this  selection  of  directions.  The 
properties  investigated  included  the  tensile  strength,  elastic  constants,  and  fracture  energies.  The 
tensile  strengths  were  evaluated  using  an  ideal  tensile  test,  which  provides  the  theoretical  tensile 
strength,  and  by  fitting  the  calculated  data  to  universal  binding  energy  relationships  (UBERs), 
which  permit  the  extrapolation  of  the  calculated  results  to  arbitrary  length  scales.  The  results 
demonstrate  the  ability  of  an  UBER  that  accounts  for  fracture  surface  relaxation  to  yield  a 
realistic  estimate  of  the  tensile  strength  of  a  10  micron  thick  sample  of  using  data  obtained 
through  calculations  on  nanoscopic  systems.  We  predict  that  C^CE  will  fracture  most  easily  in 

the  [Oil  (3/  2)(a/c)~  2]  direction,  with  a  best  estimate  for  the  tensile  strength  of  386  MPa  for 

a  10  micron  grain,  which  is  quite  consistent  with  flexural  strength  measurements  for  chromia. 
The  grain  become  considerably  stronger  at  the  nanoscale,  where  we  predict  a  tensile  strength 
along  the  same  direction  of  32.1  GPa  for  1.45  nm  crystallite.  The  results  also  provide  insight  into 
the  origin  of  the  direction  dependence  of  the  mechanical  properties  of  Cr2C>3,  with  the  differences 
in  the  behavior  along  different  directions  being  related  to  the  number  of  Cr-0  bonds  supporting 
the  applied  tensile  load.  This  work  was  published  in  N.  J.  Mosey  and  E.  A.  Carter,  “Ab  initio 
LDA  +  U  Prediction  of  the  Tensile  Properties  of  Chromia  across  Multiple  Length  Scales,”  J. 
Mech.  Phys.  Solids  57,  287  (2009). 

C.  Shear  Response  of  Chromia 

We  then  applied  the  ab  initio  LDA  +  U  (U-J=7.7  eV)  to  calculate  shear  strengths  of  chromia. 
First  we  needed  to  develop  and  test  a  renormalization  scheme  describing  the  behavior  of 
materials  subjected  to  shear  up  to  the  point  at  which  slip  occurs.  The  model  we  devised  is 
analogous  to  the  tensile-load  UBER  and  leads  to  a  size-dependent  definition  of  the  shear 
strength,  xc.  The  quantities  that  enter  into  this  definition  can  be  derived  directly  from 
calculations  on  nanoscale  samples,  and  thus  this  model  provides  a  means  of  predicting  tc  at 
mesoscopic  length  scales  from  data  obtained  through  atomistic  calculations.  The  model  was  used 
to  evaluate  tc  of  chromia  (CUO,)  for  three  relevant  slip  systems.  The  ab  initio  LDA  +  U 
calculations  clearly  demonstrate  the  ability  of  the  model  to  account  for  the  size-dependence  of  xc 
of  this  material.  For  instance,  the  ab  initio  LDA  +  U  calculations  yield  values  of  tc  ranging  from 
19.4  to  29.4  GPa,  depending  upon  the  slip  direction,  slip  plane,  and  level  of  theory,  which  is 
consistent  with  theoretical  estimates  of  the  ideal  shear  strength  of  the  material.  Meanwhile, 
values  of  tc  between  189  to  342  MPa  are  obtained  when  the  quantum  mechanics  data  are  used  in 
conjunction  with  the  proposed  scheme  to  predict  xc  for  a  grain  that  is  10  microns  thick.  The 
direction  and  plane  most  likely  to  slip  first  is  predicted  to  be  [123  1  ](0 1  12).  Although  we  were 
unable  to  find  experimental  data  against  which  to  compare  the  results,  ceramics  generally  have 
shear  strengths  in  the  range  of  a  few  100  MPa.  This  indicates  that  the  rescaling  and  extrapolation 
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scheme  employed  here  is  suitable  for  estimating  mesoscopic  shear  strengths.  The  development 
of  this  new  universal  scaling  law  for  slip  and  its  application  to  chromia  was  published  in  N  J. 
Mosey  and  E.  A.  Carter,  Acta  Materialia,  57,  2933  (2009). 

D.  The  PI  was  asked  to  write  a  Perspective  article  in  Science  for  a  special  section  on  theoretical 
chemistry  which  appeared  in  August  2008.  The  article  concerns  challenges  and  frontiers  in 
quantum  mechanics  based  materials  modeling.  The  paper  offers  a  candid  evaluation  of  methods 
ranging  from  the  most  accurate  (quantum  chemistry  and  quantum  Monte  Carlo)  to  the  least 
expensive  yet  often  accurate  enough  (density  functional  theory),  including  also  theories  for 
electronic  excited  states  in  condensed  matter.  The  paper  offers  recommendations  as  to  which 
methods  should  be  used  for  various  materials  classes  and  phenomena.  It  points  out  pros  and 
cons  of  each  method,  which  methods  should  not  be  used  or  at  least  when  one  may  get  either  the 
wrong  answer  or  the  right  answer  for  the  wrong  reason.  The  article  ends  by  discussing  a  variety 
of  materials  and  phenomena  that  still  pose  great  challenges  in  terms  of  accurate  and  efficient 
quantum  based  materials  modeling,  primarily  due  to  the  difficulty  of  integration  of  theories  for 
complex  materials.  (E.  A.  Carter,  “Challenges  in  Modeling  Materials  Properties  without 
Experimental  Input,  ”  Science,  321,  800  (2008)) 

E.  Mechanical  Properties  of  Iron  Oxides 

Stress-corrosion  cracking  of  steel  is  a  ubiquitous  phenomenon  in  which  steel  is  chemically 
corroded,  followed  by  mechanical  failure  induced  by  applied  loads.  The  products  of  steel 
corrosion  in  moist  air  at  high  temperature  (>570  °C)  are  Fe203,  FeO,  and  Fe304.  The  oxidation 
states  and  coordination  of  Fe  in  these  three  types  of  iron  oxides  are  different:  Fe2  and 
hexacoordinate  for  FeO,  Fe3  and  hexaco ordinate  for  Fe203,  mixed  valence  of  Fe2  and  Fe3  and 
a  combination  of  hexacoordinate  and  tetracoordinate  for  Fe304.  The  crystal  lattices  also  differ,  as 
shown  in  Figure  1:  FeO  is  a  rock-salt  crystal,  Fe203  is  a  corundum-type  crystal,  and  Fe304  has  a 
spinel  structure.  Both  electronic  and  geometric  structural  differences  will  contribute  to  trends  in 
mechanical  properties  of  the  three  oxides. 


Figure  1.  (a)  FeO,  (b)  Fe203,  and  (c)  Fe304. 


We  carried  out  a  detailed  study  of  the  tensile  and  shear  response  of  these  three  iron  oxides  using 
our  ab  initio  DFT+U  theory.  The  U-J  parameters  to  be  used  in  the  DFT+U  calculations  for  the 
iron  oxides  were  evaluated  in  our  earlier  work  (N.  J.  Mosey,  P.  Liao,  and  E.  A.  Carter, 
“Rotationally-Invariant  ab  initio  Evaluation  of  Coulomb  and  Exchange  Parameters  for  DFT  +  U 
Calculations,”  J.  Chem.  Phys.,  129,  014103  (2008)).  They  are:  3.7  eV  for  FeO  and  4.3  eV  for 
Fe203.  For  Fe304,  the  previous  unrestricted  Hartree-Fock  (UHF)  theory  scheme  used  for 
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evaluating  U-J  is  not  applicable  because  some  of  the  Fe  ions  are  of  mixed  valence  character, 
meaning  that  one  electron  is  delocalized  over  two  Fe  atoms.  Since  the  UHF  method  derives  a  U-J 
assuming  the  electrons  for  which  U-J  are  being  derived  are  localized  on  one  atom,  we  cannot 
apply  it  cases  where  the  electrons  spread  out  over  more  than  one  atom.  Instead  we  use  4.3  eV  for 
Fe3  ions  present  in  Fe3C>4  (akin  to  the  Fe3"  ions  in  Fe2C>3)  and  using  the  average  of  the  Fe2  and 
Fe3+  values  determined  for  FeO  and  Fe2C>3  (4.0  eV)  for  the  sites  that  have  mixed  +2  and  +3 
oxidation  states  in  the  crystal. 

Our  previously  derived  universal  binding  energy  relationship  for  tensile  loading  was  employed  to 
extrapolate  atomic  scale  ab  initio  DFT+U  data  to  macroscopic  sample  sizes.  The  extrapolated 
macroscopic  predictions  are  consistent  with  experimental  measurements.  The  ordering  of  tensile 
strengths  was  found  to  be  FeO  <  Fe304  <  Fe203,  which  correlates  with  increasing  ionicity  (as 
given  by  the  formal  charge  on  the  Fe  cations)  in  the  three  oxides.  The  direction  that  has  weakest 
tensile  strength  is  predicted  to  be  along  [Oil  (3/2)(a/c)22]  for  Fe203  and  along  [001]  for  both  FeO 
and  Fe304.  The  direction  dependence  of  tensile  properties  for  these  three  iron  oxides  can  be 
understood  via  a  local  bond  strain  analysis.  We  also  predicted  that  loading  FeO  along  [110]  and 
Fe304  along  [001]  or  [111]  produces  a  plastic  response  prior  to  brittle  fracture  at  high 
temperature.  This  work  was  published  in  P.  Liao  and  E.  A.  Carter,  “Ab  initio  DFT+U 
predictions  of  tensile  properties  of  iron  oxides,  ”  J.  Mater.  Chem.,  20,  6703  (2010). 

Turning  to  shear  response,  we  used  the  renormalization  model  proposed  in  our  earlier  work  on 
chromia  to  extrapolate  atomic  scale  shear  strengths  (derived  from  ab  initio  DFT+U  calculations) 
to  those  of  macroscopic  sample  sizes.  Like  the  tensile  strengths,  the  extrapolated  macroscopic 
predictions  are  consistent  with  measurements.  The  shear  strengths  of  FeO  and  Fe304  are 
predicted  to  be  similar,  -4.5  MPa,  while  the  shear  strength  of  Fe203  is  predicted  to  be  the  largest, 

-7  MPa,  for  samples  5  mm  thick.  The  preferred  slip  systems  are  predicted  to  be  {0001}  ^1010^ 

for  a-Fe203  and  |l  10}  ^1 10^  for  both  FeO  and  Fe304.  These  shear  strengths  are  significantly 

lower  than  the  corresponding  tensile  strengths,  suggesting  that  these  oxide  scales  are  likely  to 
slip  before  fracturing.  This  work  was  published  in  P.  Liao  and E.  A.  Carter,  “Ab  initio  DFT+U 
predictions  of  the  shear  response  of  iron  oxides,  ”  Acta  Materialia,  58,  5912  (2010). 

These  findings  have  implications  for  understanding  wear  of  gun  tubes,  since  shear  forces 
dominate  when  the  projectile  move  upon  gun  firing.  It  suggests  that  any  area  where  exposed 
oxide  has  formed  may  lose  the  oxide  due  to  shear  followed  by  mass  removal  (hence  wear). 
Moreover,  our  predictions  suggest  that  polycrystalline  oxides,  especially  on  the  nanoscale,  may 
be  the  most  wear-resistant  due  to  their  improved  strength. 
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